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ABSTRACT 

We present a revised metallicity distribution of dwarfs in the solar neighbourhood. 
This distribution is centred on solar metallicity. We show that previous metallicity 
distributions, selected on the basis of spectral type, are biased against stars with solar 
metallicity or higher. A selection of G-dwarf stars is inherently biased against metal 
rich stars and is not representative of the solar neighbourhood metallicity distribution. 
Using a sample selected on colour, we obtain a distribution where approximately 
half the stars in the solar neighbourhood has a metallicity higher than [Fe/H]=0. 
The percentage of mid-metal-poor stars ([Fe/H]<-0.5) is approximately 4 per cent, in 
agreement with present estimates of the thick disc. 

In order to have a metallicity distribution comparable to chemical evolution model 
predictions, we convert the star fraction to mass fraction, and show that another 
bias against metal-rich stars affects dwarf metallicity distributions, due to the colour 
(or spectral type) limits of the samples. Reconsidering the corrections due to the 
increasing thickness of the stellar disc with age, we show that the Simple Closed-Box 
model with no instantaneous recycling approximation gives a reasonable fit to the 
observed distribution. Comparisons with the age-metallicity relation and abundance 
ratios suggest that the Simple Closed-Box model may be a viable model of the chemical 
evolution of the Galaxy at solar radius. 
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1 INTRODUCTION 

The aim of this paper is to present and discuss a new metal- 
licity distribution of the stellar material found in the solar 
neighbourhood. In recent years, several such metallicity dis- 
tributions, constructed from various sa mples of solar neigh- 
bourhood dwarfs have been published (|Rocha-Pinto & Ma- 



ciel, 19961; |Wyse fc Gilmore, 1994 |Flynn fc Morell, 1997} 
Favata et al., 1997[). These studies have repeatedly pointed 



to a deficit of metal-poor stars relative to the simplistic but 
insightful Simple Closed-Box model (hereafter SCB model). 
However, since the early seventies, models have successfully 
fit the observed metallicity distribution by_alleviating one 
or another assumption of the SCB model D. Such fits are 
now routinely obtained with present day models o f chemical 
evolution of the Milky Way (Chiappini et al., 1997; Prantzos 
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The assumptions of the SCB model are that the solar neigh- 
bourhood is considered as closed box, with no matter flowing in 
or out, the gas is initially free of metals, the initial mass function 
is constant, and the interstellar medium is well mixed at all times. 



& Silk, 1998) using a pletora of solutions, the most widely 
accepted being the infall model. 

In view of the apparent easiness of these models to fit this 
constraint of chemical evolution, one may wonder why the 
deficit of metal-poor stars is regularly addressed, and why 
we intend to do so. One reason is that the completion of 
Stromgren surveys by Olsen (see Olsen (1983) and there- 
after) has considerably extended the available data. More- 
over, the coincidental publication of the Hipparcos Cata- 
logue now allows a clean definition of a complete sample of 
solar neighbourhood dwarfs. 

A second reason is that the role of the thick disc is still 
unclear and difficult to evaluate : is the thick disc the first 
epoch of the disc formation ? Is it cogenetic to the Galaxy ? 
And how should it be taken into account when considering 
the G dwarf problem ? In the present study, we take the view 
that the thick disc is an integral part of the disc and that 
it should be considered when discussing the G-dwarf prob- 
lem. This raises the question of the exact contribution of the 
thick disc at [Fe/H]<-0.4. Wyse & Gilmore (1995) have con- 
tended that the important metallicity tail below this limit 
(>20%) is likely due to contamination by the thin disc. This 
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problem is then connected to another. If the thin disc contri- 
bution at [Fe/H]<-0.4 is as important as envisaged by Wyse 
& Gilmore (1995), it leads to the conclusion that the SFR 
must have been very efhcient at early times in the disc. This 
conclusion is reinforced by the apparent rapid rise of metal- 



licity in the age-metallicity relation (Scully et al., 1997 



leaving even less time to produce the material at [Fe/H]<- 
0.4. As a matter of fact, it is a feature common to most 
models with infall that they use a decreasin g SFR, in accor- 
dance with the infall decay rate ( Tosi, 1996 ). In contrast, all 
recent determinations of the SFR history of the disc point to 



a co nstant, or even slight ly increasing SFR ( Haywood et al 



1997 



Binney et al., 2000). We show here that these apparent 



contradictions are mostly an effect of pre-Hipparcos data, 
and that the metallicity tail at [Fe/H]<-0.4 is drastically 
reduced with Hipparcos parallaxes. 

Third, while chemical evolution models have increased 
considerably in complexity over the last decade, the sophis- 
ticated solutions envisaged to solve the 'G-dwarf problem' 
(such as infall, variable initial mass function, etc), are still 
essentially adhoc, and direct observational evidence favoring 
one or another alternative have remain extremely elusive. 
Infall models have focused most of the efforts of galactic 
chemical evolution studies in recent years, while relatively 
little attention has been devoted to other solutions. Since 
the G-dwarf problem remains the main (indirect) support 
for long time-scale gas accretion by the Milky Way disc, it 
is important to check the failure (or success) of no-infall 
models. Also, when considering scenarios where the thick 
disc is envisaged as a genuine galactic population (as is the 
case here), it must be born in mind that no-infall models 
may come as a more natural alternative. The existence of 
the thick disc, if considered as cogenetic, implies that a sub- 
stantial stell ar disc must have been in place at early times 
(Wyse, 2000 |astro-ph/0012270| ). As remarked by Wyse, this 
suggestion could find some support or invalidation in a more 
general context, from studi es of high redshift galactic discs 
(Brinchmann fc Ellis, 2000[). 



Finally, another incentive for this work has been the 
(naive) realisation by the author that while models of chem- 
ical evolution predict distributions of stellar mass as a func- 
tion of metallicity, previous studies of the dwarf metallicity 
distributions have only been able to give statistics of stars 
as a function of metallicity. The accurate positioning of the 
stars in the HR diagram allowed by Hipparcos parallaxes 
permits the conversion from (colour, magnitude) to masses. 
It will be demonstrated that, while the conversion itself has 
a relatively minor effect, the colour limits of the sample 
bias the distribution by underestimating the contribution 
of solar-metallicity and metal-rich stars. 

The remainder of the paper is divided into 5 parts. In 
the following section, we present the calibrations used to 
estimate metallicity from Geneva and Stromgren photom- 
etry. After selecting a basic sample from considerations of 
completeness and available photometry, we select long-lived 
stars and discuss the biases introduced by our selection pro- 
cedure. A first estimate of a corrected metallicity distribu- 
tion of star frequency distribution is given. This distribu- 
tion is then compared to previous works in section 3, where 
we also discuss the problem of biases in other samples. We 
evaluate the proportion of X-ray emitter stars in the sam- 
ple. Since X-ray coronal activity is usually considered as 



a tracer of young stars, we discuss the significance of the 
high percentage of young star candidates in the sample. In 
section we convert our dwarf metallicity distribution to 
a mass metallicity distribution. A proper correction for the 
mass bias and a final metallicity distribution are given for 
both iron and oxygen abundances. In section n, we briefly 
review the observed parameters of the disc that enter the 
SCB model, and following Occam's Razor, we explore how 
the Simple model fits other local constraints of chemical evo- 
lution. We conclude in the last section. 



2 THE OBSERVED DWARF METALLICITY 
DISTRIBUTION 

2.1 Description of the sample 

The solar neighbourhood, as sampled by the Hipparcos Cat- 
alogue, is considered to be essentially complete to A'f i,=8.5 
for stars with tt >40 mas (Jahreiss & Wielen, 1997), and 
contains 959 entries within these limits. At M„=8.5, main 
sequence stars have B — V colour between 1.3 and 1.5. 
Present day available photometric metallicity calibrations 
are limited to bluer colour. The Geneva photometric sys- 
tem is probably the best suited system in that rega rd, since 
the c alibration is available down to B2 — 1^1=0.65 (Grenon, 
corresponding to B 



197S 



V^=1.05 approximately. Accord- 
ing to the available data, the main sequence of the galactic 
globular cluster M92 ([Fe/H]~-2) passes thro ugh the point 
B-V=1.00, M„=8.0 (IStetson & Harris, 198^). This means 



that the limit at M„=8.5 is unlikely to produce any signifi- 
cant bias against low metallicity stars in the sample. 

Within these limits, the sample contains 681 entries. If 
we further select stars with Mi, >3.5 and B — V >0.25 and 
which are not flagued as either G, O, S, or V suspected bina- 
ries, we find 475 entries, of which 82 are flagued "C" compo- 
nent solutions in the main catalogue. For 27 of these 82 sys- 
tems, no Geneva photometry is a vailable from the Genera l 
Catalogue of Photometric Data (Mermilliod et al., 1997). 
Those "C" component solutions for which photometry is 
available are usually not separated or simply not detected 
as binaries in the GCPD data base (except for 3 systems, 
for which photometry for the primary is available). We de- 
cided therefore to exclude entries with the flag "C" from 
further consideration, keeping in mind that we are however 
excluding a minimum of 82 stars with B — V <1.05 from 
our sample. We are then left with 393 stars for which iron 
abundance estimates are necessary. 



2.2 Photometric metallicity calibration 

When dealing with the metallicity distribution of stars in the 
solar neighbourhood, one is interested by long-lived stars 
down the main sequence, in order to avoid bias favoring 
young - and possibly more metal-rich - objets. However, 
metallicity measurements and calibration of M stars are still 
in infancy, and one is limited to stars bluer than M spectral 
type. There is no one single metallicity indicator for all the 
stars in the samples analysed in Sec. 3 and 4 and we rely on 
Stromgren or Geneva photometric metallicity, and for a few 
stars on spectroscopic data. We discuss herebelow the pho- 
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Figure 1. Pho tometric determ ination of metallicity using the 
calibration by (Grenon, 197a[), versus spectroscopic measure- 



ments for dwarfs with 0.40< B2-V1 < 0.65. 



Figure 2. Strdm,gren metallicities versus spectroscopic vaetallic- 
itie s. Photometric metalliciti es are derived from, the calibration 
by (Schuster & Nissen, 198i). 



tometric calibrations we use, by reference to spectroscopic 
measurements. 



2.2.1 Geneva photometric calibration for dwarfs : Grenon 
(1978) 

A metallicity calibration of Geneva colour indices exists for 
main sequence stars down to B2 — Vi=0.65 {B — 1/=1.05) 
(Grenon, 1975). Towards hotter main sequence stars, this 
calibration is valid up to B2 — Vi=0.4. The metallicity is 
related to the index (5i256, defined as the difference between 
the U — B Geneva colour index of the star and the Hyades se- 
quence at the B2 — V1 colour of the star (see Grenon (1978)). 
This calibration is : 



[Fe/H]Gen=2.96-F2.04/(5i256-0.72). 



Numerous spectroscopic metallicities have been pub- 
lished for G and K dwarfs since 1978, so that we can check 
how this calibration behaves compared with recent spectro- 
scopic metallicity measurements. We gathered a sample of 
spectroscopic metallicities from the literature, selecting stars 
t hat have Geneva photometry and 0.40<B2-Vi < 0.65 from 
(Edvardsson et al., 1993; Feltzing fc Gustafsson, 1995 ; [Plynn 



jFlynn 



fc MoreU, 1997| ; |Axer et al., 1995| ; [Castro et al., 1997[ ). We 
added a sample of Hyades stars, with individual spectro- 
scopic [Fe/H], as listed in Ferryman et al. (1997). The total 
sample amounts to 76 stars. 

Fig. nl shows the correlation between the photometric 
and spectroscopic iron abundance in this sample. Disper- 
sion increases for low metallicity stars, however the photo- 
metric metallicity still nicely correlates with spectroscopic 
metallicities, with no apparent large deviations. A least 



square fit to the points of Fig. nl yields a regression curve 
of [Fe/H]Gene^a = 1.002*[Fe/H]spectr-o-H0.054, indicating that 
the photometric iron abundance may overestimate the spec- 
troscopic iron abundance by an amount of 0.05 dex. 



2.2.2 Stromgren metallicity scale 

Most previous studies of the G-dwarf problem are based on 
Stromgren photometry, it is therefore interesting to see how 
Stromgren metallicity calibration compares with respect to 
the Geneva photometric metallicity and to the spectroscopic 
scale. We apply the calibration given by Schuster & Nissen 
(1989) : 



\Fe/H] = 1.052 - 73.21mi + 280.9(6 - y)mi + 
333.95(6 - y)ml - 595.5(6 - yfrm + (5.486 - 

41.62mi - 7.963(6 - y){log{mi - C3)) (1) 

when 0.22< b-y <0.375 



[Fe/H] = -2.0965 + 22.45mi - 53.8mi - 62.04mi6 -y + 

145.5mf (6 -y) + (85.1mi - 13.8ci - 137.2m?) * a (2) 
when 0.375< b-y <0.59. 

The indix c3 is defined as c3 = 0.6322 - 3.58(6 - y) + 
5.20{b-yf 

Fig. ti shows the spectroscopic vs Stromgren photomet- 
ric determination of [Fe/H] for an enlarged sample of 110 
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2.3 Age and mass biases 

One of the predicted quantities of cliemical evolution models 
of tlie galactic disc is the fraction of stellar mass as a function 
of metallicity. Ideally, a distribution also function of the age 
would be desirable, but since measurements of the age of 
(most) long-lived stars is impractical, such a distribution is 
still out of reach. However, because we observe metallicity 
to be a rough function of age, an unbiased sample must 
not favor any particular age, and in particular should not 
include short lived stars, which are only representative of 
recent chemical evolution. 

A second bias, that has been neglected in previous de- 
terminations of the metallicity distribution, is introduced in 
observed samples by the fact that, due to selection criteria 
(e.g, colour limits or spectral types), the metallicity is not 
sampled over the same width of mass interval. The selection 
of long-lived stars and the colour limit imply that metal-rich 
stars are selected on a smaller mass interval than metal-poor 
stars. This would be of no consequence if chemical evolution 
models take this effect into account, but this is never done 
in practice. We review these 2 biases in turn. 



Figure 3. Stars with Stromgren photometric iron abundance es- 
timates (hashed histogram) and Geneva photometric abundance 
estimates (thick line histogram). 



stars, with the Stromgren metallicity calculated using the 
Schuster & Nissen (1989) calibration. A fit applied to this 
sample yields [Fe/}i]stromgren=0.865[Fe/R]spectro-0.052, in- 
dicating a problem in the calibration, a result similar to 
Alonso et al. (1996). The correction to this calibration sug- 
gested by these authors is 0.85 [Fe/H]-0.04, which is concor- 
dant with our fit. 

Of the 393 stars of our sample, 324 have Geneva pho- 
tometry, 69 stars have not. Out of these 324 stars, 177 have 
0.4 < B2-V1 < 0.65, a nd their metall icitv can be estimated 
with the calibration of ( Grenon, 1978 ). Among the 147 stars 
with B2-V1 out of this interval, 9 stars have B2-V1 > 0.65, 
and are not considered any further. We are left with 138 
stars plus 69 stars which have no Geneva photometry. Out 
of these 207 objects, 166 have a Stromgren photometry, and 
the corrected calibration above can be applied. There are 
41 stars for which no Stromgren nor Geneva photometry is 
available, and 7 of these have a metallicity in the catalogue 
of Cayrel de Strobel et al. (1997). Two stars with a position 
above the main sequence that suggest binarity have been re- 
moved. Our final sample has 348 stars with estimated iron 
abundance. 

Fig. ra shows the photometric [Fe/H] distributions of the 
two subsamples (Geneva and Stromgren photometric metal- 
licities), and illustrates that they cover approximately the 
same [Fe/H] range, with a tendency for the 'Geneva' sample 
to be more metal-rich. The colour criteria for the selection 
of the stars in each sample explains this difference, the ob- 
jects in the 'Geneva' sample are redder, and it is expected 
to favor metal-rich stars. 



2. 3. 1 Correcting for age bias 

The selection of long-lived dwarfs is necessary to avoid bias 
favoring young stars in the sample, or recent galactic evo- 
lution. We are looking for those stars in the sample which 
would still be below the turn-off after a duration equal to the 
disc age. If the (thick and thin) disc age is 12 Gyr, this im- 
ply that the stars we are interested in have a main sequence 
mass lower than the turn-off mass of the 12 Gyr isochrone 
at the star metallicity, which we note M^q . In principle, be- 
cause the stars in the sample may have any age between 
and 12 Gyr, they must be located to the right of the evolu- 
tionary track with mass MJ^q , and to the left of the 12 Gyr 
isochrone. 

In practice, these requirements are impossible to apply 
because (1) we don't know the age of the disc to a satisfac- 
tory accuracy, (2) although excellent, Hipparcos parallaxes 
are not sufficiently accurate to locate meaningfully the ob- 
served star between the isochrone and evolutionary track 
(3) theory of stellar evolution and the various transforma- 
tions necessary to compare the observed star with isochrones 
and evolutionary tracks are still too imprecise. As a con- 
sequence, we have applied a simpler procedure which con- 
sisted in selecting stars to the right or below the evolutionary 
track with mass M^q . The turn-ofl' mass is obtained from a 
12 Gyr isochrone, which is taken here as the age of the galac- 
tic (thin and thick) disc. Lebreton (2000) has demonstrated 
that standard stellar evolution calculation shows differences 
of about 100-200 K compared with the best available data 
for deficient and mildly deficient stars. According to Lebre- 
ton (2000), these differences can be eliminated if (1) non- 
ETL corrections are applied to the metallicity data and (2) 
microscopic diffusion of heavy elements is accounted for in 
stellar models, which results in effective temperature shift of 
about 100-200 K. Before selecting the stars, we artificially 
applied an equivalent shift in metallicity to the observed 
sample, and in temperature to the stellar models. 

The selection of long-lived dwarfs, applied to the 348 
initial stars, yielded 218 objects. 
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Figure 4. Description of the sample of long-lived dwarfs (a) Iron abundance distribution. The upper histogram (thick line) is the sample 
before the selection of long-lived dwarfs. The (thin line) histogram counts the stars left after having selected long-lived stars, (b) HR 
diagram of selected stars, (c) For a given star in the sample, the abscissa gives the B — V of the solar isochrone at the magnitude of 
the star, minus the B — V colour of the star. 



Fig. Ma shows the iron abundance distribution result- 
ing from this selection. The initial set of 318 stars is cen- 
tred on the solar metallicity, and the distribution of long- 
lived stars (218 objects) is almost symetrically centred on 
[Fe/H]=-0.05-0.0, a slight and expected shift. Although the 
selection preferentially removes metal-rich stars, note that 
the effect is relatively minor. The metallicity histogram of 



long-lived dwarfs has 38 per cent objects with [Fe/H]>0.0, 
43 per cent for the initial sample. The Fig. Wb illustrates 
the position of the stars in the Hertzsprung-Russel (HR) di- 
agram, and Fig. Wc shows more distinctly the variation of 
metallicity along the main sequence width. 
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Figure 5. Mass interval between our two limits defining our sam- 
ple on the main sequence (Turn-off and B — V =1.05), for differ- 
ent metallicities. Crosses give the interval in solar masses. The 
diamonds give the fraction of mass interval relative to the mass 
interval at [Fe/H]=-0.7. Correcting factor has been applied after 
this. 



Figure 6. KineTnatically unweighted metallicity distributions of 
long-lived dwarfs. The hatched histogram is our brute metallicity 
distribution after selection of 218 long lived dwarfs. The thick 
line histogram is the distribution corrected for the mass bias, as 
described in the text. 



2.3.2 Approximate correction of the mass bias 

As such, the resuhing metaUicity distribution of 
Fig. W(a) should not be compared with chemical evolu- 
tion models, because the colour limits of the sample (at 
B — 1/=1.05 and turn-off) introduce a bias in the mass in- 
terval over which stars of different metallicity are sampled. 
Chemical evolution models work with stellar mass percent- 
age at a given metallicity, without any reference to how this 
mass is displayed over 'real' stars. However, at B — 1^=1.05, 
stars of metallicity [Fe/H]=-0.70 have masses of 0.60 Mq, 
whereas at the same colour, solar metallicity stars have 0.63 
M0, and at [Fe/H]=-|-0.4 stars have masses of 0.79 M©. 

Fig. y shows how the width of the mass interval 
(Mass(Turn-off)-Mass(i3 — V=1.05)) varies a s a function of 
met allic ity (crosses 
1994|)7}Vi 



in stellar isochrones by (Bertelli et al. 
i.s can be seen. 



stars with metaUicity [Fe/H]=-0.70 
are sampled over a mass interval which is roughly two times 
larger than those with [Fe/H]=-|-0.4. A correction can be 



appliedl to tne number distribution ot t ig. 5 a, wUicn is ot 
course deptiudMiL ou Lha amouul of aUi'H craaLyd &i a givaii 
mass for a given metallicity. Our correction is a linear in- 
terpolation to the fraction of the mass interval at different 
metallicity relative to the mass interval at [Fe/H]=-0.7. This 
correction is plotted on Fig. ^ (diamonds). 



2.4 The metallicity distribution of long-lived 
dwarfs 

The effect of the mass correction on the metallicity distribu- 
tion of long-lived dwarfs is shown on Fig. H. While strictly 
speaking the resulting distribution is still not comparable to 



chemical evolution models (it is not a distribution of mass 
fraction), it is an unbiased distribution of star numbers as 
a function of metallicity. We postpone the derivation of our 
final metallicity distribution of stellar mass to section W. 

The main characteristic of our distribution is that it 
is centred on solar metallicity, while previous studies have 
found distributions centred between [Fe/H]=-0.3 and -0.1. 
In order to illustrate that our result is not an effect of a 
flawed metallicity calibration, we plot on Fig. W the HR 
diagram of our stars divided into objects with metallicity 
[Fe/H]<-|-0.14 and [Fe/H]>-|-0.14, which limit is the meta, l- 
licity of the Hyades according to (Ferryman et al., 1998). 
Single stars within 10 pc of the H yades cluster centre as 
selected by (Ferryman et al., 1998) are also plotted. The 
figure clearly shows that according to their position in the 
HR diagram these stars are undoubtedly metal-rich stars. 

There are two other reasons that suggest we do not 
overestimate the global shift of the distribution to higher 
metallicities. T he first one comes from a study by ( Thevenin 
& Idiart, 19991), which shows that due to non-LTE effects. 



the spectroscopic abundance scale should be corrected. This 
correction is unimportant for solar metallicities, but is of 
the order of 0.05 dex at [Fe/H]=- 0.5 . The second reaso n 
is studied by (Gimenez et al., 1991), (Morale et al., 1996), 



( [Favataet al., 1997| ), ( |Rocha-Pinto fc Maciel, 1998[ ), showing 
that chromospheric activity affects the photometric indice 
mi on the Stromgren scale, in the sense that young active 
stars m ay appear metal-deficient. We discuss this point in 
section 3.3 



Since previous studies on the G dwarf metallicity have 
found a relative consensus that the distribution peaks be- 
tween -0.3<[Fe/H]<-0.15, we think it is highly desirable to 
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spend some time looking for the origin of the differences in 
the metalUcity distributions. 



3 COMPARISON WITH PREVIOUS 
METALLICITY DISTRIBUTIONS 

3.1 General comments 

Different biases may affect the metallicity distribution due 
to selection criteria of the observed stars. Two such biases 
have been discussed and corrected in the previous section. 
In the literature, only the first one has been corrected (the 
age bias), although it is the least important one. An even 
more problematic bias occurs prior to the two former biases 
mentioned, in samples selected on the basis of spectral type. 
A selection on spectral type is sometimes imposed by the 
availability of photometric surveys (such as Stromgren sur- 
veys of Olsen), but seems to be perpetuated for rather his- 
torical reasons (the 'G-dwarf ' metallicity distribution) than 
real limitations in the available data. Unfortunately, spectral 
type criteria introduce biases in the metallicity distribution 
which are almost impossible to correct^ 

We have demonstrated in section 2.3.2 how the colour 
limit of the observed samples introduces a non-negligible 
bias against metal-rich stars. In order to be able to correct 
this bias, one must know the colour limits of the sample. This 
is impossible to obtain once the selection has been made 
from spectral type, due to the intrinsic high dispersion of 
colours at a given spectral type. A brief inspection of the 
Third e dition of the Catalogue of Nearby Stars (hereafter 
CNS3, (Gliese & Jahreifi, 1991)) illustrates this point: for 
example G5 V stars have colours extending from B — 1^=0.58 
to 1.1. and GO V stars from 0.50 to 0.69. The intricacy of 
colours between spectral types forbids a clean definition of 
the sample mass limits. 

The metallicity bias introduced is even more severe if 
one does not include suffici ently late spe ct ral types. Thi s 
point has been discussed in (Grenon, 197^ ), (Grenon, 1987). 
According to ( Grenon, 1990D , turn-off of the oldest solar 
metallicity stars is B — ^=0.68, which implies that samples 
selected on spectral types earlier than G2 V or G5 V would 
miss most of them, not to speak of metal-rich stars. We show 
below that an unbiased sample must comprise spectral type 
as late as K. 

In the last five years numerous studies of the local 



metallicity distrib u tion of dwarfs have been publis h ed ( Wyse 
fe Gilmore, 1995|; |Rocha-Pinto & Maciel, 1996|; |Flynn & 



Morell, 1997| ; |Favata et al., 1997| [Rocha-Pinto fc Maciel, 



199S). All are based on trigonome tric distances from the 



CNS 3 or CNS2. The first two and (Rocha-Pinto fc Maciel, 



199S) used Stromgren photometry, while Flynn fc Morel 



designed their own metallicity indicator based on the Bi 
Geneva and the R — I Cousins colour indices. Finally (Fa- 



vata et al., 1997) relied on their own spectroscopic determi- 



nations. As a general rule, it can be said that the metal- 
licity distributions of all these authors present the same 
characteristics. Briefly : a broad maximum is found be- 
tween [Fe/H] =— 0.3 and —0.1, 20-30 percent of the stars 
have [Fe/H]>0.0, 10 to 20 per cent between -1.0<[Fe/H]<- 
0.5. Very few or no stars have [Fe/H]<-1 (however Flynn 
& Morel find 7 stars with [Fe/H]<-1.0) or [Fe/H]>0.3. The 
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Figure 7. The sample of long-lived stars separated into 2 parts. 
Plot (a) shows the Hyades sequence (open symbols) and the 40 
stars in the sample t hat have IFe/Hl>+0.1 i. which is the metal- 
licity of the Hyades (Ferryman et al., 199i), plot (b) shows stars 
with [Fe/HI<-h0.14- . The Hyades stars have been selected by (Fer- 
ryman et al., 199S) as single stars within 10 pc of the cluster 



centre. The sun is located at (0.62, 4-81). These figures clearly 
illustrate the existence of a large proportion of metal-rich stars 
in the sample. 



© 2001 RAS, MNRAS 000, 1-17 



8 M. Haywood 



exception to this rule is (Favata et al., 1997), with a distri- 



bution that peaks at a somewhat higher metallicity than the 
others with the consequence that 41 per cent of their stars 
have a metallicity higher than solar, and 7 per cent between 
-1.0<[Fe/H]<-0.5. 

In order to understand the origin of the differences and 
similarities b etween our metallicity distribution and those 



obtained by (Wyse fc Gilmore, 1995; Rocha-Pinto & Maciel 



199(:; Favata et al., 1997), we examine in more detail the 



samples used by these authors. 



3.2 Comparison with other studies 



l \Wyse & Gilmore, 1993^ ) 



Out of 90 stars in the Wyse & Gilmore (1995) sample given 
as nearer than 30 pc in the CNS3 catalogue, 28 are in fact 
at a larger distance according to Hipparcos parallaxes. The 
great majority of these were selected on the basis of trigono- 
metric parallax in the CNS3 corresponding to a distance 
between 20 and 30 pc. All had large errors on parallaxes. 

An interesting point is the presence of metal-rich stars 
in our sample which are totally absent from the sample by 
Wyse & Gilmore. This difference is an illustration of the bias 
introduced by spectral type. The sample of Wyse & Gilmore 
(1995 ) resulted fr om the cross-identification of the CNS3 
and (Olsen, 1983) catalogue of Stromgren photometry for 
F/GO spectral types. Wyse & Gilmore (1995) then applied 
a 2"'' selection by excluding stars outside the interval 0.4< 
B — V <0.9. A consequence of their spectral type selection is 
that 10 stars in our sample with [Fe/H]> 0.2 and 0.4< B-V 
<0.9 are absent from the Wyse & Gilmore (1995) sample. 
Nine of these objects are in the CNS3, but have G5V-K0V 
spe ctral types, and are therefore absent from the catalogue 
of (Olsen, 1983). Even more surprising is the fact that for 39 
stars with [Fe/H]>0 and B — V <0.9 in our sample, none has 
a (Simbad) spectral type earlier than G5. Therefore, 39 stars 
selected as long-lived objects and [Fe/H]>0 in our sample 
are absent from the sample of Wyse & Gilmore (1995). 

Finally, note that due to the accurate Hipparcos paral- 
laxes, the overdensity of metal-poor stars found by Wyse & 
Gilmore (1995) is washed out. Fifteen stars (~ 6 per cent) 
have [M/H] <-0.5 in our sample, 17 (19 per cent) in Wyse 
& Gilmore (1995). Six stars only are in common. Nine stars 
in the sample of Wyse & Gilmore (1995) have a parallax 
smaller than 40 mas and have not been included in our sam- 
ple. Two stars have a parallax larger than 40 mas but have 
a multiplicity flag in the Hipparcos catalogue. 



(Rocha-Pinto & Maciel, 199t) 



(Rocha-Pinto & Maciel, 1996) have selected a sample sim- 



ilar to that of Wyse & Gilmore (1995), except that they 



■f They provide a clear illustration of a Lutz-Kelker effect on 
trigonometric parallax. We note in passing that 18 of these stars 
also have a spectroscopic or photometric parallax in the CNS3, 
which in most cases is much nearer to the Hipparcos parallax (to 
within 10 per cent for most). 



extended their selection to all G dwarfs in the Gliese cata- 
logue. Thus, they are less sensitive to bias against metal- 
rich stars (but see below): approximately 20 per cent of 
their sample has [Fe/H]>0. They find very few stars with 
[Fe/H]>0.2 (approximately 2 per cent), and their distribu- 
tion peaks between [Fe/H]=-0.2 and -0.3. Several factors 
contribute to the importance of this peak. First note that 
47 per cent of the stars are illegitimate in the sample of 
( [Rocha-Pinto fc Maciel, 1996 ) , because of parallaxes smaller 
t han 40 mas according to Hipparcos. Second, the use of 
( pchuster fc Nissen, 1989 ) metallicity calibration underes- 
timates the metallicity of the stars. Finally, the limitation 
to spectral types G is undoubtedly a source of bias, (even 
though it is less important than for the previous study). 
Fig. H(a) shows their metallicity distribution, after cleaning 



from objects with -khv 



< 40 mas. We also removed 



those objects which are subgiants and giants according to 
their position in the HR diagram. The result is not very 
different from their initial distribution (see their Fig. 2). 

On the same plot, we show our sample, from which 
we have rem oved all K-type stars, in orde r to mimic the 
selection of (Rocha-Pinto & Maciel, 1996). The effect of 



this simple modification is to produce a new distribution 
whi ch peaks at -0.3<[Fe/H]<-0.2 , in agreement with that 
of (Rocha-Pinto & Maciel, 1996). The bias is clear : by 



simply removing K-type stars, the sample is shifted from 
[Fe/H]^0.0 to -0.3<[Fe/H]<-0.2 



{Favata et al., 1997) 



The distribution obtained by Favata et al. is interesting be- 
cause (1) it is based on spectroscopic measurements and (2) 
its maximum is shifted towards solar metallicity stars, as 
in our sample. They have measured the iron abundance for 
91 stars among the 1979 edition of the catalogue of Nearby 
stars (Gliese & Jahreiss, 1979), limited to stars with n > 
0.045". The distance definition of the sample suffers from 
problems mentioned for previous samples, with 40 per cent 
of these 91 stars being outside the 22.5 pc sphere. However, 
although the sample is obviously incomplete in distance, it 
is probably not biased by spectral type because the objects 
were selected at random in a given colour interval. The re- 
sult on Fig. B(b) is a distribution very similar to that of our 
sample. 



Conclusions 

Examination of the first two examples shows that they cer- 
tainly suffer from biases introduced by a selection on spec- 
tral types. Although the attention in the literature has been 
mainly focused on biases affecting low metallicity stars, we 
argue that at least as much important biases have truncated 
the metal-rich ([Fe/H]>0) side of the metallicity distribu- 
tion. All three samples are plagued by incompleteness of 
the colour intervals and underestimated distances from the 
Catalogue of Nearby Sta rs. As a conclusion, we shoul d say 
that the two studies by (Rocha-Pinto & Maciel, 1996) and 
Wyse & Gilmore (1995) have provided G dwarf metallicity 



distributions consistent with each other. However, G dwarfs 
are not representative of the metallicity of the disc, and are 
biased against metal-rich stars. On the contrary, when the 
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Figure 8. Plot (a) illustrates how a bias against metal-rich stars 
is introduced by selecting G-type stars. The plot shows our sam- 
ple of long-lived dwarfs, with all K-type stars remov ed, to mimic 
the metallicity bias (in thick line). The sample of (Rocha-Pintc 



& Maciel, 1991), which contains only G-type stars (cleaned using 
Hipparcos parallaxes) is also shown. The lower plot shows our 
distribution (all sta rs, corrected as des cribed in section 2) and 
the distribution by (Favata et ah, 199''i), which shows no appar- 



ent bias. Histograms have been normalised to contain the same 
number of stars. 



sam ples are not limite d to G spectral types, as is the case 



for (Favata et al., 1997), the metallicity distribution contains 
about 40 per cent of atara with [Fo/H]> 0.0. 



3.3 



A glimpse on the age distribution of the 
sample 



In models describing the chemical evolution of the galactic 
disc with infall, the infall rate fixes the pace at which the 
galactic disc is built, through the SFR. This is well illus- 
trated in models reviewed by Tosi (1996). Because the infall 



rate is usually taken as a decreasing function of time, the 
SFR history follows. Thus, since infall is so widely used to 
explain the lack of metal-poor dwarfs observed in the so- 
lar neighbourhood, a desirable test of these models would 
be the age distribution of dwarfs that make the metallicity 
distribution. Unfortunately, as mentioned previously, most 
of the stars in the sample have a main sequence evolution 
that is too close to the ZAMS to have an age determined 
with some accuracy in the HR diagram. However, an inter- 
esting information is available for the young dwarfs because 
they can be detected as X-ray emitters due to coronal ac- 
tivity. The cross-identification betwee n the ROSAT survey 
and the Catalogue of Nearby Stars (Huensch et al., 199!:) 
can be used to quantify the percentage of potentially young 
stars in our sample. 

There are 90 stars (40 per cent) in our sample in 
the list of objects resulting from this cross-identification. 
Fig. H(a) shows the X-ray flux of these stars as a func- 
tion oi B — V colour. The flgure shows that our sample 
may contain slightly over 90 stars with X-ray luminosity 
brighter than log X/Li,oi=-5.5, which could be absent from 
the sample of Huensch et al. (1999), due to the incomplete- 
ness to the X-ray data for the reddest stars. In order to 
get a rough estimate the corresponding ages, we utilise the 
log Lx/Lf,o!-logR'Hif relation of Sterzik & Schmitt (1997) 
and the logRVx-logt relation of Soderblom et al. (1991). 
According to these relations, log Lx/Li,oi=-5.5 corresponds 
to an age of 2 Gyr. Most our stars having X-ray emission at 
the level of log Lx/Lftoi >-5.5 could therefore be considered 
younger than 2 Gyr. 

Fig. Hb) illustrates what the metallicity distribution 
of these stars is. It is centred on [Fe/H]=0.0, with an 
unexpected but important contribution of 11 stars with 
[Fe/H]<-0.3. Some of these stars have a spectroscopic 
iron abundance which confirms the photometric iron abun- 
dance. For instance, HIP 88622 has [Fe/H]p feoto=-0-506, 
and [Fe/H]s peeiro=-0.47, dated 15 Gyr in (Edvardsson 
et al., 1993| ), and HIP 15510 has [Fe/H ]„hoto=-0.35, and 
[Fe/H]3pectro=-0.48 (Pasquini et al., 1994). The X-ray emis- 
sion of HIP 88622 may appear somewhat puzzling because 
of its metallicity, and because HIP 88622 has no significant 
chromospheric emission according to Pasquini et al. (1994). 
While these cases may be exceptions, it has been suggested 
more generally that activity may affect photometric indices 
and that presumably young stars may be measured as de- 
f icient objects if abundance is measured from photometry 



(iGimenez et al., 199l|; Morale et al., 199e| : 



Favata et al., 1997 



Rocha-Pinto & Maciel, 1998). This effect may explain the 
shape of the histogram of Fig.bl but its importance is uncer- 
tain, and we do not try to correct it. While some of the stars 
may have been detected as X-ray emitters in the sample for 
reasons not related to their age, it is clearly demonstrated 



that X-ray emission is a tracer of young stars ( Guillout et al 



199S). In this respect, the fact that around 40 per cent of 
the sample is composed of X-ray emitters that have ages less 
than 2 Gyr suggests a few comments. 

In section H, we calculate a model metallicity distribu- 
tion assuming a constant SFR. We argue that although a 
constant SFR is not phenomenologically directly related to 
the gas content, it is compatible with available determina- 
tion of the SFR history in the Milky Way. If a constant 
star formation rate has dominated the history of the galac- 
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Table 1. Percentage of stellar mass as a function of metallicity 
in the solar neighbourhood per 0.1 dex bin 




0.0 
photometric 



Figure 9. (a) X-ray flux as given in (Huensch et al., 199i} as a 
function of B — V colour. The horizontal line shows approximate 
flux leve l for 2 Gyr stars accord ing to the logLx-logR'uK rela- 
tion o£j^St£rzit£_&_^chmitt, 199', ) and the logR'uK-logt relation 
199]^ 



of (Soderblom et al., 



The other line delineates the in- 



comp leteness of the X-ray data, according to ( Sterzik & Schmitt, 



199',). The 90 stars may slightly underevaluate the number of 



X-ray emitters in our sample, because of incompleteness of the 
X-ray data B — V >0.8. (b) Photometric metallicity distribution 
of stars in (a). 



tic disc, then of the order of 15-25 per cent of the stars are 
expected to have an age less than 2 Gyr (for a 8-12 Gyr thin 
disc). However, a distance limited sample of the solar neigh- 
bourhood is biased against old stars because of the secular 
heating of the disc, observed in the increase of the vertical 
velocity dispersion with age. Using the correction factors in- 
troduced in section 5.2.2, (factor of 5 decrease in the surface 



density for stars between 8 and 4 Gyr, 2.2 for stars between 
4 and 1 Gyr, normalised to 1 for stars with age less than 
1 Gyr. The thin disc age is taken to be 8 Gyr; adopting 
10 Gyr, this estimate changes to 40 per cent), we obtain of 
the order of 45 per cent of the local stellar material to have 
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an age less than 2 Gyr. This is compatible with 40 percent 
of the stars having X-ray emission. 



4 FINAL DISTRIBUTION 

4.1 [Fe/H] distribution 

As already mentioned, a distribution suitable for compar- 
ison with presently available chemical evolution models is 
the percentage of mass found at different metallicities. To 
obtain such distribution, we e stimate the mass of each star 
from the adequate isochrone (Bertelli et al., 1994) - for un- 
evolved stars - or evolutionary track - for those stars for 
which evolution is significant. The masses of unevolved stars 
have been estimated by matching the observed and theo- 
retical absolute magnitudes on the isochrone. We have not 
taken the colour into account because of the uncertainty of 
t he stellar model temperatures and the conversion to colour 
(Lebreton, 2000). For this reason, it is expected that for 
some stars that fall below their isochrone, the mass may 
correspond to a colour on the isochrone that is beyond the 
colour limit at B-V=1.05 shown on Fig. [Oja. Stars near the 
turn-ofT all fall within the limit, because the fit to the stellar 
models includes the age as a supplementary parameter. The 
resulting fraction of stellar mass as a function of metallicity 
is plotted on Fig. [lO|b. 

In order to correct for the bias introduced by the colour 
limit and the selection of long-lived objects, we use the fol- 
lowing procedure. We assume a given initial mass function, 
the same at all metallicities over the mass interval of in- 
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Figure 10. (a) : Mass-raetallicity distribution for dwarfs in the 
sample. The 2 lines show the selections due to the colour limit 
at B — V =1.05 (left) and turn-off (right). The gray areas shows 
the limits within which a correction has been applied, (b) : the 
final metallicity distribution of stellar material in the solar neigh- 
bourhood. The result is illustrated adopting 3 different IMF slopes 
(0.05,0.7,1.2), which shows negligible effect (the 3 curves). The 
gray area shows the uncertainty in the distribution due to poisson 



N stars in the sample contribute to a given metallicity in- 
terval A [Fe/H], we can evaluate the total mass yN stars 
would represent, weighted by the adopted IMF. The result- 
ing distribution, together with its errors is given in Table 

4.2 [O/H] distribution 

In order to have a distribution comparable with the predic- 
tions of the Simple Closed-Box model with instantaneous 
recycling approximation or other types of models which 
use the instantaneous approximation, the iron distribution 
is usually converted to an oxygen abundance distribution. 
Since Pagel (1989), this is achieved by noting that oxy- 
gen, being mostly synthesized in type II supernovae, whose 
recycling time is negligible, is better suited to compare 
with the Simple model. Clegg et al. (1981) have shown 
that [O/H]=0.52[Fe/H]-f0.03 for -l<[Fe/H]<-K0.4, which 
has subsequently been approximated to [O/H]=0.5[Fe/H] 
in studies dealing with the G-dwarf metallicity distribu- 
tion. Applying this relation to the [Fe/H] distribution much 
strengthens the narrowness of the distribution and the dis- 
agreement with the SCB model, and it deserves more atten- 
tion. There is now a large debate on the exact behavior of 
the [O/Fe] ratio as a function of [Fe/H] for [Fe/H]<-1.0. We 
are interested here in the metal-richer part at [Fe/H]>-1.0. 
When measured on the dataset by Edvardsson et al. 
(1993), a linear regression yields 



[O/H] =0.64*[Fe/J/] - 0.041 



(3) 



while the study by (Carretta, Gratton, & SnedenCar- 



retta et al., 2000) gives the relation 
[O/H] =Q.hl*[Fe/H] - 0.064 



(4) 



and the study by ( Chen et al., 2000 ) gives the relation 

[O/H] =Q.70*[Fe/H] + 0.07 (5) 

The difference of 0.11 dex offset between the relations 
ra and H is due to different temperature scales adopted in 
the two studies, as mentioned by Chen et al. (2000). In the 
metallicity range of interest, the second relation gives sim- 
ilar results to the relation of Clegg et al. (1981), implying 
a narrow [O/H] distribution, but with a shift of 0.1 dex. 
On the contrary, the relation of Chen et al. (2000) seriously 
reduces the effect of converting from [Fe/H] to [O/H]. Look- 
ing at Fig [yl, it seems that the conversion to [O/H] comes 
as an additional uncertainty in an already long list of ap- 
proximations whose effect has not been seriously accounted 
for. There are too few indications that can help us decid- 
ing which conversion to apply, hence, we prefer to keep on 
working with the [Fe/H] distribution. 



terest, from 0.60 to O.95M0. This IMF is used to extrap- 
olate the amount of stellar mass outside the mass interval 
where completeness is achieved. The number of stars as a 
function of [Fe/H] within the mass interval where complete- 
ness is achieved is used as normalisation. The plot of fig. 
|lC| b shows the resulting metallicity distribution for different 
IMF slopes. The influence of the IMF slope is negligible, due 
to the small mass interval involved. We have also evaluated 
the contribution of the poisson noise to the distribution. If 



5 THE SIMPLE CLOSED-BOX MODEL 

Two kinds of assumptions define the Simple Closed-Box 
model. The first kind assumes that complex physical pro- 
cesses can reasonably be represented by simple analytical 
laws. They reflect our (poor) knowledge of these processes 
but are convenient for the derivation of the model results. 
These assumptions include : (1) an IMF constant with time, 
(2) no exchange of matter (3) the interstellar medium is 
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Figure 11. Distribution of stellar mass wi th lO/HJ, adopting 
3 diff erent IFe/HJ to [0/HI relations from i Edvardsson et al. 



199: 



(a). I Carretta, Gratton, & SnedenCarretta et al, 200C) 



(bJTiChen et al, BOOCjJc] 



initially free of metals (4) the interstellar medium is homo- 
geneous at all times. Assumptions of the second kind are in- 
troduced purely for analytical tractability. This is the case 
for the instantaneous recycling approximation (IRA). We 
are tempted to classify assumption (4) in the second cat- 
egory, because it is possible to conceive sub-SCB models 
which would evolve independantly at slightly different rates 
and that would preserve individually the characteristics of 
the SCB model. This would give rise to some amount of 
inhomogeneity in the metallicity distribution without cor- 
rupting the Simple model as the main evolutionary path for 
chemical evolution. 



5.1 The 'G-dwarf problem' 

semantics 



A little bit of 



The 'G-dwarf problem' is usually referred to as the lack 
of metal-poor stars in the vicinity of the sun, relative to 
the SCB model. However, the 'G-dwarf problem' seems to 
cover two different problems. In Section pi we have demon- 
strated that there is a real 'G-dwarf problem in the sense 
that a sample of G-dwarf stars is inherently biased in regard 
to metallicity, and cannot be compared directly to chemi- 
cal evolution models. Perhaps the 'G-dwarf problem' should 
more appropriately only refer to this specific bias. 

The second aspect of the so-called G-dwarf problem 
(which arguably should better be quoted as the SCB model 
problem), refers to the long-lasting difference between the 
amount of metal-poor material generated by the SCB model, 
and the observed metallicity distribution. 



5.2 Model input parameters 

The derivation of the metallicity distribution of the stellar 
material in the case of the SCB model with IRA can be made 
analyti cally, due to the assumptio ns of the model, see for ex- 
ample (Binney & Tremaine, 1987). In the case of zero initial 
metallicity of the gas in the disc, the metallicity evolves as 



Z{t) = -pln[Sgas{t)/Sgas{0)] 



(6) 



where S(0) is the initial gas surface density, and S(t) is 
the gas surface density at time t. 

In the above equation, the yield p modulates the po- 
sition of the peak of the metallicity distribution, and S(t) 
is an implicit combination of the SFR consumption of gas 
and the rejection rate of the stars. It is usually chosen so 
that the final ratio Sgas(t)/S5as(0) meets the (loose) obser- 
vational constraints of the solar neighbourhood. When the 
instantaneous approximation is not assumed, the gas frac- 
tion evolves as the complex result of stellar ejecta, lifetimes, 
and gas definitively locked up in low mass stars and stellar 
remnants. The resulting metallicity distribution is a function 
of the yield and the gas fraction of the disc, while the com- 
parison with the data also depends on the disc thickening. 
We review each of these three points in turn. 



5.2.1 The yield 

The yield of a generation of stars is the mass of metals syn- 
thetized and ejected per unit mass of material locked for 
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a sufficiently long time in stars and stellar remmants. It 
is dependant on nucleosynthesis, mass loss rate and stellar 
ejecta, but also on the IMF. In the case of the SCB model 
with IRA, the yield is chosen so that the equation above is 
compatible with the present fraction of gas in the disc and 
abundance of the interstellar medium. The abundance of the 



interstellar medium is usually evaluated to be solar (Binney 



& Meriifield, 199S) (this parameter is however very uncer- 
tain, and not necessarily representative of the evolution on 
the last Gyr. ). The gas fraction is situated between 0.1-0.3, 
depending on the adopted gas density (6-llM0.pc~^) and 
total density (4O-5OM0.pc~^). With these values and equa- 
tion pi the yield is evaluated to be within the range 0.0086- 
0.025, which is a fairly large range, and is of little help in 
defining the paramaters of the SCB model with IRA. 

How does this estimate compare with the yield calcu- 
lated from the local IMF and stellar yields ? Such yields 
have been calculated using a model with the following char- 
acteristics : The gas is ejected from stars assuming the same 
initial to final mass relation as that used in Scully et al. 
(1997), which comes from Iben & Tutukov (1984) : 




0.5 1.0 1.5 

slope X between 0.1-1.0 solar mass 



2.0 



m<6.8M0 
m>6.8M0 



Mfl=0.11xm+0.45 Mq 
Mfl=1.5M0. 



The stellar yields come from Portinari & Chiosi (1999) 
for massive stars and van den Hoek & Groenewegen (1997) 
for intermediate mass stars. 

Although the IMF of solar neighbourhood stars is still 
highly uncertain, it is now admitted that it has at least 
two different regimes, for low (roughly lower than one solar 
mass), and high masses. The IMF at low masses is essential, 
but in many cases in the literature the value of the index 
adopted in chemical evolution models seems unrealistic. For 
instance, Portinari & Chiosi (1999) adopt x=1.35 for M< 
2 M©, while others Gratton et al. (2000) use Scalo (1986), 
which is also too steep at low masses. Other works utilise 
a single index for calculating the yield, such as Pilyugin & 
Edmunds (1996). Recent measurements of the IMF from 



Figure 12. Yield weighted by the IMF as a function of the IMF 
slope for low mass stars (< IMq) and high mass stars (>1Mq). 
The yield is calculated assuming stellar yields at solar metallicity, 
with specifications mentioned in the text. 



the thin disc and the thick disc. Assuming the thick disc is 
responsible for 2 per cent of the local stellar mas s, and has 
a scale height of 1400 pc (Reid & Majewski, 1993), the thin 



disc with 325 pc, then the relative amount of thick disc is of 
the order of 8-9 per cent. The total surface density of visible 
stars (i.e no stellar remnants) is of the order of 27 M0.pc~^. 
Taking into account the gas surface density and allowing 
for some stellar remnant, the total surface density is of the 
order of 40-50 M0 .pc~^. Assuming disc char acteristics as 
those proposed by (Haywood et al., 1997) and (Robin et al. 



the solir neighbourhood luminosity function have proved 199^)^ one finds that the thick disc represents 15 per cent of 



that the IMl'' index at low m asses is much sha llower than 
the Salpeter value, at x=0.05 (Reid et al., 1996). At higher 



masses, the IMF is much more uncertain (see Scalo (1998), 
Haywood et al. (1997b) for a review). We consider that val- 
ues between x=l and 2 are reasonable, and values between 
2 and 3 are not unrealistic. 

Fig. ha shows that for values of the high mass IMF index 
less than « 1.9, the yield is larger than 0.01. For Salpeter 
IMF index, values of the yield as high as 0.05 are possible. 

Assuming that the present ratio Sgas(t)/St,as(0) is 
around 0.1-0.3 and that the abundance of the interstellar 
medium is of the order 0.02-0.03, this value is in the range 
of possible yields. 



5.2.2 Surface and volume densities 

An important issue when discussing the distribution of 
metallicities is the surface densities of the stellar and 
gas components in the sola r neighbourhood. According 
to (Jahreiss & Wielen, 1997), the local stellar density is 
3.9xlO~'^M0.pc~'^. The projection to distance outside the 
galactic plane can be made assuming exponential density for 



the local stellar surface density. 

A correlated problem is the correction one has to ap- 
ply in order to convert volume densities to surface densities, 
and vice versa. This is a difficult problem, because the scale 
heights of the (thin and thick) discs are a function of age, 
which is not a quantity that can be derived for stars in the 
sample. Various attempts have been made, see in particular 
Sommer-Larsen (1991) and Wyse & Gilmore (1995). Wyse & 
Gilmore (1995) estimated a correction function of the metal- 
licity. One problem with their solution is that they use the 
same correction for stars in the interval [Fe/H] = [-0.3,-1-0.3]. 
However, even though this may seem a rather narrow range, 
there is a significant age trend over this interval (see the 
age- metallicity relation below). It implies that such stars 
may have quite different ages, hence quite different correc- 
tion should be applied. We choose to 'correct' the predicted 
distribution instead (that is, we convert model surface den- 
sity to volume density), since metallicity is an explicit func- 
tion of age in the model. This solution has the advantage 
of being consistent with the fact that the vertical velocity 
dispersion is also a known function of age. 

Quantitative estimates of the corrections as a function 
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of age have been derived using the oscillation period given 
in Wyse & Gilmore (1995). We combine this relative oscil- 
lation period with the age-vertical velocity dispersion rela- 
tion given in Gomez et al. (2001). The combination of these 
two factors gives a relative correction of the order of 13 for 
'thick disc' stars, and 5 for the oldest disc stars. Note that 
these values are much larger than that applied by Wyse & 
Gilmore (1995), and partly explains why we succeed in giv- 
ing a reasonable fit to the data. Note also that these values 
are compatible with scale height variations within the (thick 
and thin) disc. 



5.3 Local constraints on the Simple Closed-Box 
model with no IRA 

Since the previous section shows the uncertainty introduced 
by converting [Fe/H] to [O/H], we preferentially work with 
the iron distribution. This means that we have to drop the 
IRA and calculate numerically the metallicity distribution. 



5.3.1 The model 

The main characteristics of the model are given in Table H 
The stellar yields are used as described in section |5.2| for 
massive an d intermediate m ass stars. The stellar lifetimes 
come from (Pols et al., 1997), and are dependent on metal- 



Table 2. Characteristics of the SCB model of Fig.l5. 



licity. The onset of type la supernovae happens when the 
metallicity reaches [Fe/H]=-1.0. The yields for the differ- 
ent species produced by type la SN are from Nomoto et 
al. (1984) (iron ^ 0.6 Me/event) and the rate of SNIa is 
assumed to be proportional to the number of SNII, with 
SNII/SNIa=8.5. Note that we don't use a Schmidt law type 
SFR. While there are some evidences that the SFR of mas- 
sive stars may be proportional to some power of the gas den- 
sity, there are few evidences that this can be extrapolated 
to low and intermediate mass stars. Since determinations to 
date are compatible with a constant SFR history for the disc 
of the Milky Way, we use this simple prescriptionlSI. 



5.3.2 Metallicity distribution 

The delayed ejection of important quantities of gas from 
long-lived stars could enhanced the dilution of metals in 
the interstellar medium, and permit the production of more 
stars at intermediate (e.g solar) abundance. In our tests how- 
ever, this effect proved to be very minor, for the reason that 
the release of gas is spread on long time-scale by the impor- 
tant variation of stellar lifetime with abundance and mass; 
therefore, there is no sudden release of gas, and the dilution 
of metals in the interstellar medium is smoothed over long 
time scales. 

The only marked effect of the SCB model with no IRA 
is the more rapid metallicity increase at the onset of SNIa 
at [Fe/H]=-1. This translates in the metallicity distribution 
into a shift of the distribution to higher metallicities, as 



3 It is expected, and perhaps has been demonstrated that to some 
level, the SFR has not remained constant, on time scale <l-2Gyrs. 
This is unimportant for the point considered here. What we mean 
is that there is no demonstration that the SFR has decreased on 
increased systematically over 10-12Gyr. 



4OM0.pc-2 

Fe/H=0. 

Constant 3.5M0.pc~2.Gyr"^ 

x=0.05 m< 1M0, x=1.7 m> IMq 

8MQ.pc-2 

14 Gyr 



Surface density 

Initial metallicity 

SFR 

IMF 

Final gas density 

Age of the model 



can be seen at [Fe/H]«-1 (Figjl3^). The SCB with no IRA 
doesn't seems to generate more solar metallicity stars. Note 
that this somewhat contradicts the findings of Scully et al. 
(1997). However, they use very different prescriptions (stel- 
lar lifetimes independent of metallicity, IMF index=1.7 over 
the whole mass range, which they take to be 0.4< m/M© < 
100). 

The total duration of the model is 14 Gyr. In the model, 
we identify three phases corresponding (not necessarily in a 
univocal correspondance) more or less to the galactic stellar 
populations. The first phase is defined by [Fe/H]<-1.0 dex. 
We do not try to ascribe a particular stellar population (halo 
or thick disc) to this phase, and the scale height correction 
attributed to this metallicity range is arbitrary. In other 
words, we don't include this part of the distribution in our 
discussion of the SCB model. This conservative position is 
justified to some extent by (1) the fact that the character- 
istics of the metal weak thick disc (or fiattened halo ?) are 
still essentially unknown (2) the limited volume of our sam- 
ple is not well suited for studying an intrinsically rare pop- 
ulation. In the model, the duration of this phase is 3 Gyr. 
The thick disc phase starts at [Fe/II]=-1.0, and may last 2 
to 3 Gyr, ending when the metallicity reaches -0.58 or -0.45. 
Then the thin disc phase may be identified with the remain- 
ing evolution, with the metallicity rising from -0.58 or -0.48 
to [Fe/H]=-|-0.17. The disc as a whole is characterised by 
having around 11 Gyr, and a metallicity that evolves from 
[Fe/H]=-1 toJFe/H]=-h0.17. 

Figure [13(b) shows the metallicity distribution of the 
SCB model with no IRA after the disc thickening has been 
accounted for. A gaussian noise of 0.15 dex has been added 
to account for metallicity dispersion at all ages, plus 0.1 dex 
for simulating individual measurement error. The figure 
shows a quite reasonable agreement between the model and 
the observations. 



5.3.3 The age-metallicity relation 

On Fig. [141 a), we display the age-metallicity rel ation (AMR) 
of the mode l, together with the one derived by (Rocha-Pintc 



et al., 2000). The plot shows a reasonable agreement between 



the two over a large age range. Note that the two points 
cor responding to very old st ars are upper limits according 
to dRocha-Pinto et al., 26oo| ). 

The age-metallicity relation h as been the subject of nu - 
merous studies since the work of (Edvardsson et al., 1993). 



Because the cosmic dispersion in the relation is viewed as 
an important constraint for chemical evolution models, it 



has f ocused a number of comments (see ( Rocha-Pinto et al. 



200C) for a review). The d ispersion in the AMR measured 
by ( [Edvardsson et al., 1993 ) is abou t 0.2 dex, larger than th e 
one given by more recent studies. (Edvardsson et al., 199J) 



cautioned however that their sample has been chosen for 



© 2001 RAS, MNRAS 000, 1-17 



A revision of the solar neighbourhood metallicity distribution 15 



0.15 



0.10 - 



0.05 - 



0.00 




tinction made (section 5.2) in the scale height corrections 
for stars with metaUicity varying between -0.5 and +0.2 is 
justified. 

There are three main processes by which the dispersion 
may be explained: inh omogeneity in the interstellar medium 
([Malinic et al., 199S), diffusion of stellar orbits (Grenon 



1987 




Figure 13. Plot (a) shows two Closed-Box models, one with IRA 
and yield p=0.02 (thin line) and no IRA (thick line), as described 
in the text. Plot (b) shows the SCB model (no IRA), converted 
to volume density as described in section 5.2.2, to which we have 
added a dispersion of 0.1 dex to account for metallicity measure- 
ment dispersion in the observed distribution. 



representativeness of all metallicities, which doesn't mean 
completeness at all metallicities. This makes their sample 
ill-suited for a correct measurement of the real dispersion. 
Although they provide corrections for the bias, their proce- 
dure is necessarily uncerta in. It i s not surprising that (G!ar- 
nett fc|Kobulnicky, 2000[) and ([Rocha-Pinto et al., 200C|y 



find a smaller dispersion, of the order of 0.10-0.15 dex. A 
consequence on the AMR is that the correlation between 
age and metallicity is ti ghter than has been thought after 
(Edvardsson et al., 199S). This view confirms that the dis- 



Francois fc Mattcucci, 1993|; [Wielen et al., 1996D , and 



sporadic infall ([van den Hoek fc de Jong, 1997). While it is 



almost certain that the first two play a role, even if minor, 
infall has the capability to explain large dispersion in the 
age- metallicity relation. 

Whatever the processes that lead to variations of the 
conditions of the star birth, it is worth noting that a 10 per 
cent dispersion on the IMF, SFR and SNl rate is sufficient 
to account for the observed dispersion. This is illustrated in 
Fig.Ma), where each point has been simulated as the result 
of SCB models with slightly different properties simulated 
at random from the model of table H allowing for 10 percent 
dispersion in the IMF slope (for m> IMq), SNIa rate and 
SFR. A w eighti ng according to the disc thickening described 



in section 



5.2.2 has been applied. Fig. Il4b shows the metal- 



licity histogram for these simulated points, compared to the 
metallicity distribution of Fig. pi 



5.3.4 Abundance ratios 

Our aim here is to demonstrate further how a model with 
simple prescriptions gives a satisfactory fit to complemen- 
tary local constraints such as abundance ratios. 

Fig.na displays these ratios for CNO and Si, Mg ele- 
ments. The model has obvious failures (N/Fe, Mg/Fe) and 
reasonable success (O/Fe, C/ Fe), in comparable proportion 
to more sophisticated models (Timmes et al., 1995|; [Portinari 



et al., 1998; 



Goswami & Prantzos, 2000D. Possible reasons 



for the failure o f N/Fe, Mg/Fe are discussed in ( Goswami & 



Prantzos, 200C), and we don't replicate their discussion. 



6 DISCUSSION 

A reanalysis of the solar neighbourhood metallicity distri- 
bution has brought the following results : 

(1) The metallicity distribution of long-lived dwarfs is 
centred on [Fe/H]~0, not [-0.3,-0.1]. When considering a 
sample of stars with no bias on the age, between 40 and 
50 percent of stars have metallicity higher than [Fe/H]=0. 
Most previous studies have used samples biased against stars 
with solar metallicity or higher. This result reconciles the 
age- metallicity relation with the dwarf metallicity distribu- 
tion. 

(2) The percentage of mid- metal-poor stars ([Fe/H]<- 
0.5) is in agreement with estimates of the thick disc density 
from remote star counts and is situated between 2-4 per 
cent. The thin disc is not an important contributor to stars 
with [Fe/H]<-0.5. 

(3) 40 per cent of these long-lived stars in our sample are 
X-ray emitters, which suggest that a large part of the sample 
may consist of young stars (age <2 Gyr). This supports the 
idea that the SFR is not a decreasing function of time. 

(4) We have evaluated the distribution of stellar mass 
with metallicity. Due to selection criteria of observed sam- 
ples, stars of varying metallicities are sampled on unequal 
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Figure 14. (a) The mean age-metallicity relation given by the 
model (thick curve). The dots are calculated according to the 
same model, but assuming 10 per cent dispersion on the IMF, 
SFR, and SNl rate, plus 0.1 dex 'measurement error' in [Fe/H] 
and age. This simulated 'local' distribution has been weighted ac- 
cording to the disc thickening described in section 5.2.2. The 
points with e rror bars are the observ ed age-metallicity relation 
according to (Rocha- Pinto et al., 200L). Note that the two points 



corresponding to ages 11.5 and 13.5 Gyr are only upper esti- 
mates, (b) Shows the histogram of metallicities for these points 
(thick line). The histogram of Fig. 6 is also shown (thin line). 
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Figure 15. Abundance ratios for Carbon, Nitrogen, Oxygen, Sil- 
icon and Magnesium. 
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mass intervals. We show that if this effect is left uncorrected, 
a second important bias against metal-rich stellar mass is in- 
troduced. 

(5) We argue that a reasonable fit to the [Fe/H] distri- 
bution can be obtained with a SCB model with no instanta- 
neous recycling approximation. Satisfactory fits are also ob- 
tained for the age-metallicity relation and abundance ratios. 
These results suggest that the solar neighbourhood metal- 
licity distribution contains few (or no) indications that long 
time-scale infall may have play a major role in the building 
of the galactic disc. 

In their study, Wyse & Gilmore (1995) commented that 
'even were the Simple Box model to fit some dataset, its in- 
herent implausibility means that it is more likely that sev- 
eral compensating effects had generated this agreement by 
chance'. How much unrealistic is the SCB model presented 
here ? Our model incorporates all ingredients usually found 
in chemical evolution models of the Galaxy, (i.e metallicity- 
dependent stellar yields and lifetimes, empirical IMF and 
SFR), and satisfies local constraints. The only marked fea- 
ture absent in our model is infall. While infall is prevalent 
in chemical evolution models, it is merely viewed as an ad- 
ditional free parameter to control the dilution of metals in 
the interstellar medium. In this respect, the SCB model is as 
much realistic as other models, but it is simpler. However, 
the remark of Wyse & Gilmore (1995) still holds, and the 
assumptions of the SCB models should be looked at with 
even more scrutiny. 

The SCB model is basically a description of the chem- 
ical evolution of the mean parameters of the galactic disc. 
What we seek in fitting the observed distribution presented 
here is to know if the evolutionary path followed by the 
SCB model is the correct one. This approach has a chance 
of being sound if the assumptions of the SCB models do not 
contradict observed features. We shortly review the main 
assumptions of the SCB model. 



troduced in the model while keeping with the SCB model 
as the main evolutionary path for chemical evolution. For 
example, it may be conceived that the disc as a whole has 
behaved as an ensemble of sub-SCB models with slightly 
varying properties. This would not seriously affect the pic- 
ture of the SCB model as the progressive enrichment of a 
gaseous disc initially free of metals. 

(4) Closed Box. This arguably is the most challenging 
hypothesis of the SCB model. Note that what we have as- 
sumed with this hypothesis is that the surface density of the 
model has remained constant, therefore fixing the rate of di- 
lution of metals. In sophisticated models where the Closed- 
Box hypothesis is not assumed, infall is essentially viewed as 
a way to regulate the metal enrichment, in order to provide 
a solution to fitting the observed metallicity distribution. It 
works essentially as an additional free parameter. In view 
of the work presented here, infall looks unnecessary to de- 
scribe the solar vicinity metallicity distribution. The local 
fossil signature of infall is still to be found. We note that 
the potentially most serious indication that infall processes 
may have played an important role in the building of the 
disc lies in the age-metallicity relation. 

Infall provides the only mechanism to generate old 
metal-rich (solar-like) stars at the solar radius. If it can be 
demonstrated that old stars (age>8 Gyr) with disc-like kine- 
matics and solar-like abundance exist in significant propor- 
tion, then they may represent the signature that the disc, 
or part of the (old) disc has formed through an evolutionary 
path that differs from the SCB model evolution. 



(1) Initial metallicity Fe/H=0. It was not necessary to 
adopt a non-zero initial metallicity to account for the ob- 
served metallicity distribution. If correct, this implies that 
stars at arbitrarily low metallicities with disc kinematics 
must exist. Morrison et al. (1990) found that thick disc 
stars possibly reach as low metallicities as [Fe/H]=-1.60, 
while Martin & Morrison (1998) have been able to show that 
objects with thick disc kinematics can reach [Fe/H]=-2.05. 
It remains to be demonstrated that thick disc metallicities 
can reach even lo wer values. Note that at the same time, 
( Fuhrmann, 1998| ) showed that the halo and thick disc may 
be contemporary. 

(2) Constant IMF. There is to date no clear-cut evi- 
dence for systematic variations in the IMF. The case for 
a univers al IMF has been considered by different authors 
(Kroupa, astro-ph/0009005; Gilmore (1999)) and is favored 
at low masses (<1Mq). The diagnostic for a universal IMF 
i s much less clear for masses of concern in nucleosynthesis 
(3calo, 1998). Variable IMF have been envisaged in the con- 
text of chemical evolution but has received relatively poor 
support (e.g Martinelli & Matteucci (2000) and Chiappini 
et al. (2000) for recent attempts). 

(3) Homogeneous interstellar medium. As mentioned 
above, this is a minor assumption of the SCB model, intro- 
duced for analytical tractability. Inhomogenity can be in- 
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Table 3. The list of 348 stars. Selected stars are marked with a "*" 
HIP [Fc/H] Mt, B -V HIP [Fe/H] M^, B ~V 



100925 


-0.02 


5.17 


0.72 




39342 


0.14 


5.99 


0.87 


* 


101997 


-0.60 


5.53 


0.72 


* 


40118 


-0.39 


5.14 


0.68 


* 


102264 


-0.37 


5.19 


0.67 


* 


40693 


0.03 


5.47 


0.75 


* 


103859 


0.11 


6.25 


0.97 


* 


40774 


-0.34 


6.16 


0.90 


* 


105038 


0.02 


6.87 


1.02 


* 


41926 


-0.43 


5.95 


0.78 


* 


105152 


0.15 


6.81 


0.99 


* 


42074 


0.04 


5.63 


0.79 


* 


105905 


-0.37 


6.82 


0.92 


* 


42499 


-0.21 


6.29 


0.83 


* 


106696 


-0.36 


6.29 


0.88 


* 


42808 


-0.21 


6.33 


0.92 


* 


107022 


0.02 


5.34 


0.75 


* 


43587 


0.35 


5.47 


0.87 


* 


107625 


0.09 


6.75 


0.96 


* 


43726 


0.11 


4.84 


0.66 




108156 


0.11 


6.23 


0.91 


* 


45170 


-0.45 


4.94 


0.73 


* 


109378 


0.22 


4.91 


0.77 


* 


46580 


-0.06 


6.69 


1.00 


* 


109527 


-0.03 


5.51 


0.81 


* 


46626 


0.21 


6.95 


0.99 


* 


110649 


0.03 


3.75 


0.67 




46816 


0.60 


6.50 


0.93 


* 


111888 


-0.14 


6.70 


0.94 


* 


46843 


-0.14 


5.76 


0.78 


* 


112190 


0.31 


6.44 


0.97 


* 


49366 


-0.10 


6.34 


0.89 


* 


112870 


-0.36 


6.68 


0.85 


* 


51271 


0.17 


7.01 


1.04 


* 


113421 


0.39 


4.68 


0.74 


* 


51819 


0.06 


5.68 


0.82 


* 


114416 


0.01 


7.16 


1.00 


* 


52462 


-0.09 


6.05 


0.87 


* 


115331 


0.01 


5.66 


0.80 


* 


53486 


0.09 


6.12 


0.92 


* 


115445 


-0.12 


6.35 


0.88 


* 


54426 


0.11 


6.56 


0.94 


* 


116085 


0.20 


5.63 


0.84 


* 


54704 


-0.02 


5.37 


0.76 


* 


116745 


0.17 


6.81 


0.99 


* 


55210 


-0.22 


5.57 


0.73 


* 


116763 


-0.29 


5.82 


0.80 


* 


56452 


-0.34 


6.07 


0.81 


* 


118008 


0.01 


6.51 


0.97 


* 


56829 


-0.10 


6.78 


0.98 


* 


10138 


-0.25 


5.93 


0.81 


* 


56997 


-0.15 


5.43 


0.72 


* 


10798 


-0.80 


5.83 


0.72 


* 


57443 


-0.34 


5.06 


0.66 


* 


12158 


0.12 


6.17 


0.94 


* 


57507 


-0.36 


5.23 


0.68 


* 


13402 


-0.09 


5.96 


0.86 


* 


58451 


0.30 


6.35 


0.97 


* 


13976 


0.24 


6.11 


0.93 


* 


58576 


0.33 


5.00 


0.76 


* 


14150 


0.16 


5.04 


0.70 




59280 


0.16 


5.53 


0.79 


* 


15099 


0.10 


6.09 


0.86 


* 


61291 


-0.18 


6.09 


0.84 


* 


15457 


0.06 


5.02 


0.68 




61451 


0.10 


6.19 


1.02 


* 


15510 


-0.35 


5.35 


0.71 


* 


61946 


0.05 


6.44 


0.95 


* 


16537 


-0.14 


6.19 


0.88 


* 


62229 


0.00 


6.30 


0.94 


* 


17420 


0.10 


6.36 


0.93 


* 


62523 


0.18 


5.12 


0.70 




17439 


-0.01 


5.94 


0.87 


* 


64690 


-0.06 


5.16 


0.71 


* 


18324 


-0.05 


6.21 


0.83 


* 


64924 


0.01 


5.09 


0.71 


* 


18915 


-1.60 


7.17 


0.86 


* 


65530 


0.22 


4.86 


0.74 


* 


19422 


0.14 


6.38 


0.95 


* 


66147 


0.28 


6.65 


1.03 


* 


21988 


0.15 


6.26 


0.91 


* 


66765 


-0.03 


5.97 


0.86 


* 


22122 


-0.02 


6.03 


0.88 


* 


67620 


-0.10 


4.94 


0.70 


* 


24874 


-0.11 


6.83 


1.00 


* 


67655 


-0.93 


5.98 


0.66 


* 


25421 


0.25 


6.46 


0.95 


* 


69357 


0.00 


6.11 


0.87 


* 


25544 


-0.23 


5.53 


0.75 


* 


69414 


-0.07 


5.31 


0.73 


* 


26779 


-0.01 


5.79 


0.84 


* 


69972 


0.39 


6.29 


1.02 


* 


27207 


0.15 


5.79 


0.83 


* 


70016 


0.09 


5.99 


0.87 


* 


27887 


0.07 


6.62 


0.94 


* 


70950 


0.06 


6.99 


1.03 


* 


28954 


-0.01 


5.83 


0.81 


* 


71181 


0.16 


6.61 


1.00 


* 


29271 


0.08 


5.05 


0.71 


* 


71395 


0.08 


6.39 


0.97 


* 


29525 


-0.15 


5.15 


0.66 




72312 


-0.05 


6.31 


0.89 


* 


29568 


-0.16 


5.26 


0.71 


* 


72688 


0.04 


6.67 


1.04 


* 


32010 


0.17 


6.89 


1.02 


* 


73005 


-0.22 


5.88 


0.79 


* 


33690 


0.12 


5.50 


0.79 


* 


74702 


-0.06 


5.96 


0.83 


* 


33852 


0.36 


6.44 


0.99 


* 


75253 


0.42 


6.27 


0.97 


* 


34414 


-0.15 


6.59 


0.91 


* 


75722 


0.15 


5.97 


0.87 


* 


36210 


-0.05 


4.96 


0.69 


* 


75829 


-0.11 


5.62 


0.80 


* 


38228 


-0.12 


5.21 


0.68 


* 


76375 


0.48 


5.91 


0.95 


* 


38784 


-0.20 


5.40 


0.72 


* 


77358 


0.06 


5.09 


0.71 


* 


39064 


0.25 


5.89 


0.83 


* 


77408 


-0.23 


5.77 


0.80 


* 


39157 


-0.69 


5.87 


0.72 


* 


78775 


-0.50 


5.87 


0.73 


* 
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Table 3. 


(continued) 














HIP 


[Fc/H] 


Mt, 


B-V 




HIP 


[Fc/H] 


M„ 


B-V 


78913 


-0.09 


6.74 


0.96 


* 


105184 


-0.06 


4.87 


0.80 


79190 


-0.35 


6.32 


0.86 


* 


105858 


-0.70 


4.39 


0.77 


79248 


0.43 


5.35 


0.88 


* 


107350 


-0.04 


4.64 


0.75 


79537 


-1.16 


6.84 


0.81 


* 


107649 


0.03 


4.60 


0.60 


80366 


-0.09 


6.73 


0.95 


* 


109422 


0.09 


3.58 


0.69 


81813 


-0.13 


5.63 


0.77 


* 


109821 


0.05 


4.51 


0.61 


82588 


-0.18 


5.49 


0.75 


* 


110109 


-0.27 


4.69 


0.51 * 


83389 


-0.27 


5.48 


0.73 


* 


112117 


-0.15 


4.13 


0.59 


83541 


0.30 


5.29 


0.81 


* 


113357 


0.33 


4.52 


0.62 


83990 


-0.22 


6.71 


0.89 


* 


113829 


-0.01 


4.72 


0.64 


85042 


0.04 


4.84 


0.68 




114886 


-0.16 


6.15 


0.49 * 


85235 


-0.44 


5.90 


0.76 


* 


114924 


-0.10 


4.04 


0.59 


86796 


0.35 


4.21 


0.69 




114948 


-0.14 


4.07 


0.60 


87579 


-0.19 


6.52 


0.94 


* 


115147 


-0.48 


6.05 


0.49 * 


88972 


0.07 


6.17 


0.88 


* 


116416 


-0.10 


6.05 


0.65 * 


90790 


-0.21 


6.20 


0.86 


* 


116613 


0.08 


4.76 


0.61 


91438 


-0.32 


5.29 


0.67 


* 


117712 


-0.29 


6.19 


0.58 * 


92858 


-0.17 


6.08 


0.86 


* 


118162 


0.15 


4.80 


0.67 


92919 


-0.70 


6.41 


0.91 


* 


12114 


-0.04 


6.50 


0.62 * 


93858 


0.04 


4.98 


0.70 


* 


12444 


0.02 


4.12 


0.87 


93966 


0.18 


4.48 


0.70 




12653 


0.18 


4.22 


0.56 


95319 


0.12 


5.42 


0.80 


* 


14286 


-0.22 


4.84 


0.52 * 


95447 


0.40 


4.25 


0.76 




14632 


0.29 


3.94 


0.89 


96085 


0.06 


6.26 


0.92 


* 


15131 


-0.40 


4.82 


0.85 * 


96100 


-0.14 


5.87 


0.79 


* 


15330 


-0.21 


5.11 


0.67 


96183 


0.10 


5.37 


0.75 




15371 


-0.22 


4.83 


0.98 


96901 


0.12 


4.56 


0.66 




15442 


-0.18 


5.09 


0.69 * 


97944 


0.38 


5.43 


1.02 




16852 


-0.17 


3.60 


0.92 


98130 


0.05 


7.00 


1.02 


* 


17147 


-0.82 


4.75 


0.52 * 


98677 


-0.29 


5.73 


0.71 


* 


17378 


0.12 


3.74 


0.56 


98767 


0.31 


4.74 


0.75 


* 


18267 


0.06 


5.23 


0.63 


98792 


-0.28 


6.32 


0.81 


* 


18859 


0.09 


3.96 


0.60 


98828 


0.22 


6.19 


0.92 


* 


19076 


0.08 


4.78 


0.59 


99240 


0.36 


4.62 


0.75 




19233 


-0.09 


4.54 


0.64 


99711 


0.13 


6.27 


0.94 


* 


22263 


0.12 


4.87 


0.60 


99825 


0.03 


6.01 


0.88 


* 


22449 


0.07 


3.67 


0.65 


1031 


-0.20 


5.68 


0.77 


* 


22451 


-0.16 


6.22 


0.57 * 


1292 


-0.03 


5.36 


0.75 


* 


23437 


-0.54 


5.28 


0.55 * 


1499 


0.14 


4.61 


0.67 




23693 


-0.04 


4.38 


0.92 


3093 


0.26 


5.65 


0.85 


* 


23835 


-0.08 


3.91 


0.72 


3206 


0.29 


6.18 


0.94 


* 


24786 


-0.04 


3.98 


0.52 


3535 


0.56 


6.32 


0.98 


* 


25278 


-0.04 


4.17 


0.62 


3765 


-0.04 


6.38 


0.89 


* 


26394 


0.13 


4.35 


0.64 


3850 


-0.37 


5.78 


0.77 


* 


27072 


-0.00 


3.83 


0.63 


3979 


-0.26 


5.21 


0.66 


* 


27435 


-0.21 


5.01 


0.48 * 


4148 


-0.07 


6.41 


0.94 


* 


27913 


0.04 


4.70 


0.90 


6379 


-0.10 


6.02 


0.83 


* 


28267 


-0.05 


5.16 


0.64 * 


6917 


-0.25 


5.91 


0.97 


* 


29432 


0.01 


5.03 


0.53 


7339 


-0.06 


4.92 


0.69 


* 


29800 


-0.07 


3.58 


0.66 


7576 


-0.05 


5.80 


0.80 


* 


30314 


0.02 


4.68 


0.57 


7734 


0.08 


4.98 


0.69 




30503 


0.14 


4.65 


0.54 


8102 


-0.50 


5.68 


0.73 


* 


30630 


-0.32 


5.95 


0.60 * 


8362 


0.02 


5.64 


0.80 


* 


32366 


-0.19 


3.99 


0.48 


9269 


0.20 


5.18 


0.77 


* 


32423 


-0.12 


6.81 


0.64 * 


544 


0.02 


5.41 


0.75 


* 


32439 


-0.09 


4.18 


0.59 


100017 


-0.13 


4.69 


0.60 




32480 


0.10 


4.15 


0.72 


101345 


0.12 


3.74 


0.69 




33277 


-0.13 


4.55 


0.64 


102040 


-0.16 


4.82 


0.61 




33537 


-0.32 


5.02 


0.43 * 


103389 


-0.03 


4.09 


0.51 




34017 


-0.11 


4.53 


0.61 


103458 


-0.78 


4.85 


0.59 


* 


34567 


0.10 


5.14 


0.63 


104436 


-0.42 


5.06 


0.62 


* 


35136 


-0.33 


4.41 


0.94 
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Table 3 


. (continued) 














HIP 


[Fc/H] 


M„ 


B-V 


HIP 


[Fe/H] 


M„ 


B-V 




36439 


-0.33 


3.86 


0.47 


75181 


-0.26 


4.83 


0.64 


* 


36515 


-0.16 


4.97 


0.64 * 75809 


-0.25 


4.85 


0.67 


* 


36704 


-0.14 


6.21 


0.86 * 77052 


0.12 


5.03 


0.68 




37349 


-0.14 


6.42 


0.89 * 77257 


0.11 


4.07 


0.60 




38908 


-0.23 


4.54 


0.57 


77801 


-0.49 


4.86 


0.60 


* 


40035 


-0.17 


3.77 


0.49 


78072 


-0.23 


3.62 


0.48 




40843 


-0.23 


3.84 


0.49 


79578 


0.31 


4.85 


0.65 




41484 


0.04 


4.63 


0.62 


79672 


0.21 


4.76 


0.65 




42333 


0.13 


4.87 


0.65 


80337 


0.13 


4.82 


0.62 




42438 


-0.14 


4.86 


0.62 


80686 


0.00 


4.49 


0.56 


* 


42697 


-0.11 


6.36 


0.90 * 81300 


0.00 


5.82 


0.83 


* 


43557 


-0.08 


4.66 


0.64 


81520 


-0.44 


5.33 


0.62 




43797 


0.11 


3.79 


0.48 


83006 


0.07 


6.71 


1.00 


* 


44075 


-0.88 


4.16 


0.52 


83601 


0.00 


4.45 


0.58 


* 


44897 


0.11 


4.54 


0.58 


84862 


-0.35 


4.59 


0.62 


* 


45333 


-0.05 


3.72 


0.61 


85653 


-0.53 


5.47 


0.74 


* 


45617 


0.19 


5.98 


0.99 * 85810 


0.22 


4.65 


0.64 




47080 


0.35 


5.16 


0.77 * 86736 


-0.16 


3.64 


0.47 




47592 


-0.09 


4.07 


0.53 


86974 


0.31 


3.80 


0.75 




48113 


0.23 


3.75 


0.62 


88348 


0.34 


5.31 


0.80 


* 


50075 


-0.06 


4.60 


0.59 


88622 


-0.51 


4.86 


0.61 


* 


50384 


-0.55 


4.03 


0.50 


88694 


-0.14 


4.74 


0.62 




50505 


-0.23 


5.09 


0.65 * 89474 


-0.01 


4.52 


0.64 




50921 


-0.27 


5.20 


0.66 * 89805 


0.09 


4.37 


0.58 




51248 


-0.39 


4.56 


0.61 * 93185 


-0.24 


4.95 


0.61 




51933 


-0.39 


3.76 


0.53 


96395 


-0.03 


4.81 


0.64 




52369 


-0.10 


4.94 


0.62 


96895 


0.16 


4.32 


0.64 




53721 


0.13 


4.29 


0.62 


97675 


0.05 


3.68 


0.56 




57939 


-1.66 


6.61 


0.75 * 98470 


-0.13 


4.05 


0.50 




61053 


0.03 


4.49 


0.57 


98505 


-0.19 


6.25 


0.93 


* 


61317 


-0.21 


4.63 


0.59 


98959 


-0.23 


4.83 


0.65 


* 


62207 


-0.61 


4.75 


0.56 * 99137 


0.02 


4.43 


0.53 




63366 


-0.27 


5.93 


0.77 * 99461 


-0.28 


6.41 


0.87 


* 


64394 


0.14 


4.42 


0.57 


1598 


-0.40 


4.99 


0.64 


* 


64457 


-0.11 


6.01 


0.93 * 1599 


-0.12 


4.56 


0.58 




64550 


-0.08 


4.99 


0.64 


1803 


0.19 


4.84 


0.66 




64583 


-0.29 


3.62 


0.49 


3497 


-0.24 


4.85 


0.65 


* 


64792 


0.04 


3.92 


0.58 


3909 


-0.11 


4.22 


0.51 




65515 


-0.02 


5.59 


0.80 * 5862 


-0.03 


4.08 


0.57 




67275 
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